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EXIT TIMES DENSITIES OF BESSEL PROCESS 


GRZEGORZ SERAFIN 


Abstract. We examine the density functions of the first exit 
times of the Bessel process from the intervals [0,1) and (0,1). First, 
we express them by means of the transition density function of the 
killed process. Using that relationship we provide precise estimates 
and asymptotics of the exit time densities. In particular, the re¬ 
sults hold for the first exit time of the n-dimensional Brownian 
motion from a ball. 


1. Introduction 

The Bessel process (t) with index /i G M is a well known diffnsion 
on [0, cxo). For indices of the form fi = n/2 — 1, n = 1, 2, 3,..., it may 
be described as a norm of n-dimensional Euclidean Brownian motion. 
The process is also closely related to the geometric Brownian motion 
by the Lamperti relation. Furthermore, in form of the Bessel-Brownian 
diffusion, it plays a crucial role in approach to hyperbolic Brownian 
motion BE] and is used to represent symmetric stable processes as 
well na. Finally, it is exploit to model a stock prices and has further 
applications to mathematical hnance BEllIS]- 

In the context of Bessel process, the killed process was recently of 
special interest BElElElDroilllllll]. In the paper, we examine the exit 
time density of the Bessel process starting from x > 0 and killed when 
it reaches a hxed level a > x. Due to a scaling property it is enough 
to consider a = 1. Apparently, there are in the literature noticeably 
more results concerning starting point greater than the killing level. 
Reasons of this disproportion are not clear. Bessel process killed when 
exiting [0,1) is interesting not only from probabilistic but also from 
analytical point of view, since its transition density is well known as 
a Fourier-Bessel heat kernel. One of obstacles when dealing with this 
case may be its complex nature. Technical complications appear even 
for /i = —1/2 (see [IS]), h is when the Bessel process corresponds to the 
one-dimensional Brownian motion. Our aim is to provide asymptotics 
and uniform estimates of the exit time density in whole range of index 
yii G M. Analogous results are already known in case x > a B US] HZ], 
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but are obtained by usage of completely different methods. Note that 
presented in the paper result hold, as a special case, for the first exit 
time of the Brownian motion from a ball and, up to the knowledge of 
the author, they were not known before. 

One of the difficulties to deal with is the behaviour of the Bessel 
process at zero, which depend on the value of the index p. For /i ^ 0 
the process never riches zero, and for /i ^ —1 we impose that the point 
zero is killing. In case /i G (—1,0) we may consider one of both; killing 
or reflecting condition at zero. Let be the first hitting time of the 
Bessel process at the point a ^ 0. Additionally, we introduce the first 
exit time ^ = Ta A Tf, from the interval (a, b), where 0 ^ a < x < b. 
According to the dual nature of the point zero we define for /i > — 1 
and reflecting condition at zero 

x) = (Tb e dt), a: e (0,1), f > 0. 


Furthermore, when p < 0 and zero is killing we denote 
9^(4. 1 , a) = r (r„,j e *; «('■> (r„,j) = y), 


where x G (a,fe), y G {a,b} and t > 0. In 1980, J.T. Kent yielded in 
his paper [El the following series formula 


( 1 ) 


(m) 

^1 


(x,f) = x 




x 


n=l 






Similar results may be obtained also for the function {t, x, y), which 
is presented in Section 3. Those formulae are convenient when dealing 
with large times. We focus on the other case, i.e. when t is small. 
The above-given sum is oscillating and, in that case, the cancellations 
between the terms really matter in the context of asymptotic behaviour. 

The main idea of our approach to the problem is to express the exit 
times densities in terms of the density of the killed Bessel process. This 
kind of formulae seems to be quite general, but the mentioned singu¬ 
larity of the process at zero requires a separate consideration. The 
obtained representations enable us to use some results and methods 
provided in na, where the density of the Bessel process killed when 
hitting the level 1 was estimated. One of the our key tools is an imita¬ 
tion of a reflection principle, which lets us efficiently approximate the 
density of the killed process for some range of parameters. Properties 
of Bessel functions are also intensively exploit. 

The paper is organized as follows. In Section 2 we introduce some 
notation and gather relevant facts concerning Bessel functions and the 
Bessel process. Section 3 is devoted to representations of the densities 
of the first exit times and theirs estimates. The precise asymptotic 
behaviour is provided in Section 4. 
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2. Preliminaries 

2.1. Notation. For two positive functions / and g we denote / < 
g whenever there exists a constant c > 0 such that f < eg for all 
arguments. If / ^ and ^ / we write f ^ g. Additionally, if the 
constant c depends on an a parameter, we put this parameter over the 
sign < or Furthermore, for x > 0 we use notation f{x) = Of^{x), 
which throughout the paper means that there exist positive constants 
Cl = Ci(/i) and C2 = C2(/u), such that |/(x)| < CiX for |a:| < C2. 


2.2. Bessel functions. The Bessel function of the hrst kind J^{z) is 
given by the series 

One can deduce from the dehnition the asymptotic behaviour of the 
function at zero 


( 2 ) 


lim J^{z)z ^ = 1/F (/i + 1). 


Combining this with the formula 8.451(1) in |S] we obtain 

fi 

(3) ILWI := - 

For /i > —1 the function Jfj,{z) has inhnitely many real zeros 
which have the following properties (see [18], p. 199, 506) 


)M+1/2 ■ 


(4) 

(5) 


1- JiiA 

hm = TT, 

fe—loo h 

(-IX' [2 

L+iX.O - „ O 


1 + 0. 



The modihed Bessel function of the first kind is given by I^{z) = J^{iz). 
The function is positive for z > 0 and its asymptotic behaviour is given 
by (cf. |B], formulae 8.445 and 8.451(5)) 


(7) Ii,{z) = ^ (1 + (^^)) ) ^^0. 


This implies 


{z) ^ 




( 8 ) 


(1 + ^)l‘+l/2 


z > 0. 
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Furthermore, there are provided in |T3] the following bounds for ratio 
I^{y)/I^{x) of modihed Bessel functions of the hrst kind 

4 ( 2 /) 


(9) 


( 10 ) 




< 




> 


/i > 1/2. 


where 0 < x < y. However, we will need in the sequel results of that 
kind valid for all y > —1. We prove therefore inequalities which are 
slightly weaker but hold in the whole required range of parameters. 

Lemma 2.1. For y > —1 and 0 < x < y we have 


yY+\y_^ ^ l^{y) 






^y-x ^ ^ 




-) 

yj 


\ V + 4: 




Proof. The proof is based on the equality (see j8] 8.486 1.) 


( 11 ) — Iy+2{,z) + 

Using it and the formula ( 


2(/i+l) 




z ^ C — {O}, y G 


L(») = L+2(!/) + kL±l)/, 


y 


< 


y \ 

a \ ..V—. 




e^-^Iy+2{x 


we get for p > —1 and 0 < x < y 

y+l{y) 

2{y + l) 


< I - I 


y \ M+2 


y-x 


'/i+2 


.X 


y 

2(/i + l) 


X 


X 


'm+i 


U+i 


X 


,x 


= U' 


yy^'^ ..V- 


y-x 


4(^)> 


which is equivalent to the hrst inequality in the thesis. The lower bound 
will be proved in two steps, since f[T0|) works only for y> 1/2 and the 
formula flTT]) let us enlarge the range of y only by 1. First, assume that 
y > —1/2. Then, by f[T0|) and flTT]) . we obtain 


hiy) = 4 + 2 ( 1 /) + 


2(/i + l) 


'm+1 


{y) 


> 


= - 


> - 


= - 


yJ 

-) 

yJ 

x\ 

y) 

x\ 

y) 


/i+2 


/1+2 


/i+2 


/i+3 






oV-x 


^/i+2 


(x) + 


2(/i + l) 


'm+1 


X 


4(^) 


y-x2{y + l) 


X 




.X 


y-x 

y 


4(^) 


oV-^ 


Iy{x). 
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Analogously we get for fi > 

\ /i+S 

x\ 




> 


y) 

x\ 

y) 


y-x 


, , 2(/i + 1) 

Ili+2{x) H- - - Ifj_+i{x) 


M+4 




[X 


This completes the proof. 


□ 


2.3. Bessel process. We denote by the Bessel process with 

index /i G M and starting from the point x G [0, cxo). It is a diffusion 
on [0, oo) with infinitesimal generator given by 

2fx + l d 


( 12 ) 




2 dx"^ 2x dx 
Let p^^\t,x,y) be the transition density function (with respect to 
the speed measure m{dx) = 2x'^^~^^) of the process If /x > —1 

and 0 is reflecting, we have 


(13) p^>^\t;x,y) =—{xy) 

(14) p(i^\t;0,y) = 
Applying 

(15) p^^\t;x,y)& 


x^ + y'^ 
2t 


t 


a:,!/ > 0, 


■exp (-i/V2t) , i/>0. 


{2tY+^V{p + l) 
to the above-given formulae we obtain 
-{x-yfRi 

fi > ■ 


- 1 . 


{yx +y/t' 

Consider now the Bessel process starting from a: > 0 and killed when 
hitting the level 1 . To distinguish two possible ways of behaviour of the 
process at zero, we will denote by p[^\t; x,y), p > —1, the transition 
density of the process killed only at 1 and by pQ^l{t] x,y), /x < 0, the 
transition density of the process killed at 1 and at 0 (equivalently: 
killed when exiting from the interval (0,1)). Both of those functions 
may be expressed by Hunt formula, which, in case of p[^\t; x, y), takes 
the form 

(16) p[^\t,x,y) =p^^\t,x,y) -E*[t > Tpp^^\t - Ti,l,y)], 

where f > 0, x, x/ G (0,1). An another formula formula is also known, 
which follows from general theory of compact operators: 


(17) 


(m) 

Pi 




n=l 


Estimates of p^f\t,x,y) were recently provided in |14| 




(1 -x)(l -X/)' 














EXIT TIMES DENSITIES OF BESSEL PROCESS 


6 


Notation and basic properties for the first exit times of the Bessel pro¬ 
cess were presented in Introdnction. However, some complementation 
is needed. Since (t) corresponds to the one-dimensional Brown¬ 

ian motions, we have (H 3.0.6 (a), (b) p. 212) 


(19) 


^ \x-y\ + 2k ^ / {\x-y\+2kf \ 


where x G (0,1), y G {0,1}. Moreover, the following scaling property 
holds for a > 0 and x G (0, a), i/ G {0, a}: 


(20) x) = 4 f 4’ 

\a^ aJ 

( 21 ) qo^l{t, x,y) = ^ q^^l 

’ ' \a^ a a/ 


The laws and of Bessel processes with indices /i G M and 
z/ G M, respectively, are absolutely continuous and the corresponding 
Radon-Nikodym derivative is described by 


( 22 ) 


f ds \ 

dPi'^^ V ^ V 2 4 [R(^)(s)]4 


where x > 0, and the above given formula holds Pi^^-a.s. on {Tq > t}. 
Taking u = —p, one can see that potential theory of processes with 
opposite indices are closely related. In particular, for p < 0 we have 

(23) p^o^}{t,x,y) = {xy)-‘^^^p[~^\t,x,y), 

(24) {t, X, 1) = x-^>^q[~'"\t, x). 


3. Formulae and estimates 

The discussion presented in the paper is based on the below-given 
theorem. The densities of exit times from intervals (0,1) and [0,1) 
are there represented by means of the density of the killed process. 
Note that if the considered intervals were bounded away from zero, the 
general theory of diffusion processes could be applied. 

Theorem 3.1. Let x G (0,1). For p > —1 we have 

q[^\t,x) = - 


(25) 


^pY^{t,x,l). 

dy 
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Furthermore, for jj, < 0 we have 

(26) 1) = 1), 

(i, X, 0) = lim (t, X, y) 

’ y->-o ay ’ 

(27) =-2px~^^^p^^f{t,x,t)). 

Proof. Since the density function of the hrst exit time is given as a 
derivative of the distribution function, we get 

F\t,x) = -^r(T,>t) 

(28) =-2^j jJf\t-,x,y)y'^i‘*'dy. 

By (jll), ([5]) and ([3]) we may differentiate the series ffT7)) term by term 
and the obtained function is uniformly integrable on (0,1) for t bounded 
away from zero. This allows us to switch the order of integration and 
differentiation in fl25]) . Using heat equation (cf. (IT^ i and integrating 
by parts, we get 


q{^\t,x) = 


'0 L 




^ .0 I 2/i + 1 (9 u) 


dy 


^pr{t-,x,y) + 


'0 


y^^^^-^Piit;x,y) 


(29) 


Applying 


dy 


^p[^\t-,x,y) 


y dy 

dy 


Pi {t-,x,y) 


y^^+^dy 


y 


2^+1 


y=l 

y=0 


lim 


([5]) and (l3|) to the series flTTl) we may deduce also that 
y^o -§ijPi^\t] x,y)y‘^^~^^ = 0, and this gives us fl25|) . Assume now 

that /i < 0 and 0 is killing. By fl23|) and fl2T)) . the formula fl26|) comes 
from the following 



= - X ^p[ ^\t,x,y) 

y=i 

= -x-^^ ^ (^{xyf^Po^l{t,x,y)^ 
= {i, X, 1 ) - ^Po!i i't, X, 1 ) 

= —^Pipit,x,l). 
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To prove the last part of the theorem, note that in the same way as 
(12^ we get 

(To,i e dt) = X, 0) + x, 1) 
d (/i) , . 

= - 

Thus, by (l26|) . we arrive at 

, a:, 0) = lim {t;x,y). 

’ y^o oy ’ 

Finally, fl2^ gives us 

Qo'!h't,x,0) = lim (^{xy)-^^^p[~^\t, x,y)^ 

= (^-‘^pp[~^\t,x,y) + y-^p[~''\t,x,y)^ 

= —2px~'^'^p[~^\t,x,0), 

which ends the proof. □ 


2/i+l 


y--L 


7 = 0 


An immediate corollary of the theorem is the series formulae for the 
function q^^l (f, x, y). 

Corollary 3.2. For /i < 0 and x G (0,1) we have 


(30) 

(31) 


(/i)/, IN —u \ ^ ■ d—iJ.{j—y,,nX) _j2 ^j2 

%i\^^X,\)=X T / -te , 

^—y+1 \3—y.,n} 


n=l 


g^7(f,x,0) = —^--^e ' 




\J—y+l{,j-y.,n) \ 


Proof. The formula (IHUj) is a consequence of flMj) and ([T]). By (HT} we 
have to calculate a limit of flT7|) as y tends to zero. Since the series 
is uniformly convergent for every fixed x G (0,1) and f > 0 (cf. fl3l)). 
usage of ([2]) leads to fl3T]) . □ 


Since p[ ^\t,x,l) = 0, the formula f[T8|) allows us to estimate the 
derivative {d/dy)p[^\t-, x, 1). Combining this with Theorem 13.II we may 
obtain estimates of densities q[^\t,x) and q^l{t,x,y). 


Theorem 3.3. For p > —1 we have 



^ (1 - x)(i + 

(x + t)l‘+l/2t3/2 


exp 
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For jj, < 0 we have 


1) 

9o!i 0) 


V 


V 


X ^^(1 — x) (1 + t) 

{x + t)-/^+V2 ^2 

x~^^(l — x) (1 + / 

^^-IZTT^-exp 


(^_1 A 

2t ) 

x"^ 1 .2 \ 

^ - 2^-wi^J ■ 


1 


Proof. The last assertion of the theorem follows directly from fl27|) and 
(HSl). Moreover, since fl24p holds, we need only to prove the first asser¬ 
tion. We get from the previous theorem that 


q[^\t,x,l) = 
We apply (fT8|) and, by 


dy 


■p[^\t-,x,l) 


lim 


p[^\t-,x,y) 

1 -y 


lim- - - 

y^i 1 — y 


1 

t’ 


the proof is complete. 


□ 


Let B{t) be the n-dimensional Brownian motion starting from x G 
M”. By q^{t, x) we denote the density function of the exit time of Brow¬ 
nian motion from the unit ball i?(0,1) centered at the origin. Since 
|i?(f)| is a Bessel process with index n/2 — 1 (with reflecting condition 
at 0 for n = 1), we have q"'{t,x) = ^\t,x) and, consequently. 


Corollary 3.4. For x G B{0, 1) we have 


q-{t,x) 


n 


(1 — |x|)(l -I- 

(|a;| +t)(n-l)/2^3/2 


2t 


22n/2-l,l^ 


4. Asymptotics 

In this section we provide precise asymptotic behaviour of the func¬ 
tions q[^\t,x) and q\f^l{t,x,y) for small t. In the Theorem 14.11 we 
describe this behaviour by means of simpler objects, i.e. the standard 
transition density p^^^ {t, x, 0) of the Bessel process reflected at zero, and 
the exit time density q^^^y^\t,x,y) of the one-dimensional Brownian 
motion from (a:/4,1). 

Theorem 4.1. For p > —1 we have 

(32) q[^'^ {t, x) = P^^^ {t, x,l)(^l + Of, () > 

(33) <,P{t,x)= (l + 0^(t)). X 6 (1/2,1), 
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and for jj, < 0 we have 

(34) (t, X, 0) = (t, x, 0) (1 + O/, ) , 

(35) qQ^l{t,x,0) = — -p^~^\t,x,0) 

ij 

X ^1 + 0^ ^ ^ ^ ’ 

where t < Iq for some to > 0. 


Proof. Using Strong Markov property we obtain for x G (0,1), p > —1 
and A C (1, oo) 


[i?^(t) G A] = E^ [Ti < t; R^^{t) G A] 

= E" [Ti < t; E^ - s) G , 

which can be expressed in the following way 

r rt 


I p'^^\t,x,y)m{dy) = I q[^\s,x) I p^^\t — sAiy)‘^{dy)ds. 

Ja Jo Ja 

By the Fubini theorem, we deduce 

p^^\t,x,y)= [ q[^\s,x)p^^\t - s,l,y)ds, x G (0,1), i/> 1. 

Jo 

This let us rewrite (IT^ into the form 

(36) p[^^ {t, X, y) = (t, X, y) - p^>^^ (t, x,2-y) 

- [ q[^\s,x,i)[p^^\t-s,i,y)-p^^Ht-s,i,‘J-y)]ds. 

Jo 

This equality may be understood as a kind of imitation of the reflection 
principle. The appearing integral plays then a role of error term and 
would vanish if x, y) was a probability density function of Brow¬ 

nian motion. Furthermore, by flT^ we have for n > 0 and x < 1 < y 


p(^)(M,x,2 - y) _ f y Y If,{x{2-y)/u) _ 2 (i_ 


2 - 2 / 


(37) , w 

pWi^u, X, y) 

Hence, by Lemma [2.11 we get 

y V , . P^^KuA,y) 


-y)/u 


I^ixy/u) 


- 1 < 


2-yJ ■ ^ p(A[u,l,‘J-y) 
We may write it for i/ > 1/2 as 


- 1 < 


2 - 2 / 

y 


2/i+4 


- 1 . 


1 , 2 /) 

p(A(u, 1,2-y) 


1 = 0 ( 1 - 2 /), 
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which gives us 

[ (s, X, 1) {t-s,l,y)- p^^'> {t- s,1,2-y))ds 

Jo 

= 0{l — y) / q^^\s,x,l)p^^\t — s,1,2 — y)ds 

Jo 

= 0{l-y)p^>^\t,x,2-y). 

This allows us to transform (IMl) to 


p[^\t,x,y) 

(38) = p^^^'>{t,x,y) 


l-{l + 0{l-y)) 


(39) =p’^^^\t,x,y){l + 0{l-y)) 


p^^^\t,x,2-y) 

p^iJ{t,x,y) 

p^>^\t,x,2-y) 

p(f^'){t,x,y) 


+ 0{l-y) 


Using fl37)) and Lemma l2Tl we get 


2-y 


2/i+4 


y \ • ^_ 2 il-x){l-y)/t ^ p^^\t,x,2-y) 


< 


V^^Ki,x,y) 

g-2(i-U(i-y)/i 

y 


This, together with ([38]), ll27|) and symmetry of p[^\t,x,y), proves the 
formula fl3T|) . Furthermore, the above-given inequalities give us 


1 ^-2(i-x)(i-,)/t 


P^^\t^x,y) 


y 


= 1 _ g-2(l-x)(l-?;)/t ^_ y ] I g-2(l-x)(l-y)/t 


y 


= 2 


{l-x){l-y) 


t 


+ Oij 


{l-y){l-x) 


+ 0^{1 — y). 


and similarly 

p^^\t, x,2 — y) 


1- 


P^^\t^x,y) 
{l-x){l-y) 


< 1 - 


y 


2 - 2 / 


2/1+4 




= 2 


+ On 


(1 -2/)(l - + 


+ Ofn{l — y). 
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which can be generally written as 
-y) 


P^^Kt,x,y) 

{l-x){l-y) 


= 2 


+ Or 


(1 -y){^-x) 


+ l-y 


Combining this with (15^ we arrive at 

{t, X, y) = pO) X, y) (1 + 0^{l-y)) 


X 


yi - x)(i - i/) 


+ l-y 


This eqnality implies formulae (l32|) and fl35|) . Indeed, we have 


y/'l I —y 


= 2 


1 — X 


+ 0^(1) ]pO\t,x,l). 


and 


? 0 T 0) = -2/^a; (t, 0, x) 

= — 2/1 (1 + O (1 — x)) 2^—^pO\t, 0, x) 




X I 1 + 

1 — X 


1 — X 


t 


+ t 


= —4/x " ^ pO^ (t, X, 0) ( 1 + O 


1 — X 

t 


It is showed in [T3] that there exists to > 0 such that 


(1 + 0^(t)). 


(40) p[^\t,x,y) = {l + Of,{e pi/i,i(t, x, i/), 

whenever x,y E (1/2,1) and t < to. By the absolute continuity fl22|) of 
distributions of Bessel processes with different indices we get 






= 


t<T™,/t(-V2)(t)eA; 

i?(-l/2)(i)\^+V2 / ^2_W4 

- exD--- 


X J 


ds 




Since > 1/8 on {t < we have 

< 64t, 


(a;/ 4 ,l)J 

0 ds 

lo [RF^is)? 
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and consequently 


exp 



-1/4 

2 


= 1 + 


exp 


— 1 + 0^{t). 


/i^ — 1/4 /■* ds \ 

^ Jo 



Combining this with (140]) we obtain 


p[^\'t,x,y) = (1 + 0^(t)) 

Thus, formulae fl25|) and fl26|) . scaling property and the shift-invariance 
of one-dimensional Brownian motion give us 


qf{t,x) = (1 + 0(f))^ 


d 


dy 




y=l 


(l + 0(f)) 


(/i+i/2)p^f (t, X, 1)+X, 1) 


^/i+l/2 

_ (1 +0(t)) (-1/2). . 

- ^y+l/2 ^ 4,1 

This completes the proof. 


dy^ 


□ 


We are now prepared to express the asymptotic behaviour explicitly. 


Theorem 4.2. For jj, > —1 we have 


(41) q[^\t,x,l) 

(42) q[^\t,x,l) 


1 _ T p-(l-3:)^/2i / 

1 — X / 

fM+i 2/^r(/i +1) + 


If ja < 0 , then we have 



(1 -h Ofj, (t)). 


(43) 

(44) 

gj/^(t,x,0) 


2x-^>^ 

(2f)-/^+W(-p) 




(1 + 0 , 


8(1-x)x 2/21 

(2f)-/^+2r(-p) 


1 + 0^ 



(1 + 0,(t)), 


where t < to for some to > 0. 


Remark. Asymptotic behaviour of (t, x, 1) follows from the above¬ 
given theorem and the relationship fl241) . 


Proof. To obtain fl43|) and (144)) . it is enough to apply (114)) to (l34)) and 
(135)) . respectively. Furthermore, if x/t —?■ 0, then f/(l — x) ~ f. Hence, 
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by (l32|) and ([7]), we get (142|) . What has left is to prove (HT]) . Rewriting 
(HD for y = 1 we get 


1 — X 




exp 


(1 — x)" 

2r~ 


^ exp 

k=—oo 


2k{k + 1 — x) 


2k 


^ ^g-2fc(fc+l-a:)/t _ g-2A;(fc-(l-x))/t j 


k=l 

1 — X 


exp 


(1 — x)" 

2r~ 


[S'! + S 2 ] ■ 


It is clear that S'! = 1 + 0^(e By inequality sinhr < re'", r > 0, 

we get 


2k 


0 < 5*2 = ^ ^g-2fc(fc+l-a;)/t _ g-2fc(fc-(l-a:))/t^ 


fc=l 

00 


^*^^/*sinh (2/c(l — x)/t)) 


fc=i 


^—2xlt 


< 


t 




2 -2k{k-l)/t 


k=l 


Hence, for 1/2 < x < 1 we obtain S 2 = 0^{e and we arrive at 
d-l/2)r 


= 


(1 + 0 ^(e-‘/‘)) . 


Thus, using fl33|) . scaling property and the shift-invariance of one¬ 
dimensional Brownian motion, we obtain 




V^t3/2 


^^3/2 

In addition, (l32|) and ([7]) give us for 0 < x ^ 1/2 


(l + 0^(e-"(i-"/")/*)) (l + 0^(f)) 

1/2 < X < 1. 


1 + I — 


1 ) = 


1 — X 

x/^f3/2 a;M+l/2 


1 -(- O, 


X 


(1 + (t)) ■ 


Since t/x > t, the above-given formula is equivalent to fITT]) . This ends 
the proof. □ 
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